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We propose an approach for collective enhancement of precision for remote optical lattice clocks and a way of 
generation of the Einstein-Podolsky-Rosen (EPR) state of remote clocks. In the first scenario a distributed spin 
squeezed state (SSS) of M clocks is generated by a collective optical quantum nondemolition measurement on 
clocks with parallel Bloch vectors. Surprisingly, optical losses which usually present the main limitation to SSS 
can be overcome by an optimal network design which provides close to Heisenberg scaling of the time precision 
with the number of clocks M. We provide an optimal network solution for distant clocks as well as for clocks 
positioned in close proximity of each other. In the second scenario, we employ collective dissipation to drive two 
clocks with oppositely oriented Bloch vectors into a steady state entanglement. The corresponding EPR state 
provides secret time sharing beyond the projection noise limit between the two quantum synchronized clocks 
protected from eavesdropping. An important application of the EPR entangled clock pair is remote sensing of, 
for example, gravitational effects and other disturbances to which clock synchronization is sensitive. 
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Introduction. Optical atomic clocks provide some of 
the most precise and accurate physical measurements 
to date ma. The precision of optical lattice clocks is 
presently limited by the available frequency stability of 
the best lasers EHEl , but quantum noise of uncorrelated 
atoms lures not far below. With enhanced laser stabil¬ 
ity and improved measurement protocol to reduce the 
laser noise EHm, the next frontier of precision can be 
advanced by generating entangled states of the clock 
atoms. Eor small number of atoms N, a maximally 
entangled clock operating with GHZ states can reach 
the Heisenberg limit of stability much faster than the 
best classical schemes (IS [131. Spin-squeezed states 
(SSS) m are particularly suitable for improving the 
precision of optical lattice clocks that operate on large 
N and currently hold the record for clock precision 0. 

Distant clocks connected into a spin squeezed network 
can provide a higher collective precision for all users. 

Spin squeezing (SS) and entanglement of atomic en¬ 
sembles have so far been experimentally demonstrated 
for single ensembles of spins associated with atomic 
states separated by radio- or microwave frequencies. 

This was achieved by optical quantum nondemolition 
(QND) measurement (HHsa, by mapping squeezed 
light onto an atomic ensemble 1^ , by atomic inter¬ 
actions in a Bose-Einstein condensate ED, and by en¬ 
gineered dissipation 1^ . Improvement to clock preci¬ 
sion beyond the quantum projection noise (QPN) was 
demonstrated for microwave clocks ifT^I^ . 

Networks of remote clocks offer new possibilities 
for secret time sharing, remote sensing and interferom¬ 


etry that can take advantage of the unprecedented clock 
precision. A recent proposal outlined probabilistic gen¬ 
eration of GHZ-type of entanglement by single photon 
communication between distant clocks each contain¬ 
ing a small number of Uq qubits as discrete quantum 
variables |[l2l[T3l. However, for optical lattice clocks 
containing macroscopic numbers of atoms N we en¬ 
counter continuous variables such as spin squeezed and 
EPR entangled states that can be generated determin¬ 
istically. QND probing on cyclic transitions has been 
identified as the condition for Heisenberg scaling with 
N in Ref. ll26ll and demonstrated for microwave clocks 
in Ref. (m El. However, Heisenberg scaling with 
the number of clocks M in a chain is problematic 
as SSS have low tolerance for losses, in particular to 
losses of the optical channel for optical QND. Here we 
demonstrate that an optimal design of the network of 
cavity-enhanced optical clocks allows to keep Heisen¬ 
berg scaling with M even in the presence of substantial 
channel losses. Optical lattice clocks with their long 
coherence times are ideal for generation of such states. 

In the second part of the Letter we present a sce¬ 
nario where an Einstein-Podolsky-Rosen (EPR) entan¬ 
gled state of two clocks is generated by engineered dis¬ 
sipation. Such clocks feature the ’’synchronized time” 
protected from any eavesdropper and available only for 
the participants working together. Both proposals are 
aimed at optical clocks with a macroscopic number of 
atoms. As a specific example, we show their feasibility 
for Sr clocks with realistic experimental parameters. 

A network of clocks in a collective squeezed state. 
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The ultimate limit of precision for a clock made 
of N independent atoms is defined by the the angu¬ 
lar uncertainty of a coherent spin state (CSS) of the 
ensemble (quasi-)spin CSS is a product state 

1^) = n^^^(|l)^ + \2)i) of uncorrelated atoms ori¬ 
ented in the same direction, = N/2. The other two 
projections of J have minimal equal variances allowed 
by the Heisenberg uncertainty relation: Var{Jz) = 
Var{Jy) = Jxl‘2^ = N/A. Introducing quantum cor¬ 
relations between atoms allows to reduce Var{Jz) to 
below the QPN limit. For the resulting spin squeezed 
state (SSS), under the condition that 

^ AT = Var{XA)j < 1, (D 

the atoms become entangled 1^ . and the correspond¬ 
ing signal-to-noise ratio for spectroscopy is improved 
by the inverse of ^ 03, which is the spin squeez¬ 
ing parameter for metrology. The canonical opera¬ 
tors, Xa = JzI^^Pa = Jyl\^ obey the com¬ 
mutation relation Pa] = i where J is the length 
of the mean pseudo-spin vector. The quantum noise 
limited clock precision defined as the minimal de¬ 
tectable angle of spin rotation for the clock sequence 
is yyar(J,)/J = v^. 

As shown theoretically ll^ and experimentally d 
for Cs clocks, QND probing of the clock levels using 
cyclic transitions under a high optical depth leads to 
Heisenberg scaling of precision with N. This approach 
led to the recent demonstration of SSS for Rb ground 
state 1^ . 

The QND interaction H oc XaXl leads to the 
input-output relation for photonic canonical variables 
Xl and Pl 

pout ^pin j^I^Xa. ( 2 ) 

with K = y^dr]e~^ the interaction constant, d the res¬ 
onant single pass optical depth, r] = ndr{"f/cr/A a 
parameter describing spontaneous emission caused by 
the probe, a the resonant dipole cross section, A the 
beam cross section, 7, A the natural linewidth and de¬ 
tuning of the optical transition, Udr the photon number 
in the QND probe. For QND on a cyclic transition the 
degree of spin squeezing is 

? = _n /-I 2 \ - 

e ^ (1 + K,^) 

This results in ^rnin = V^lf^opt — ^e/d oc N~^ 



FIG. 1. Clock operation and probe. The atomic level struc¬ 
ture for the optical lattice clock, (a) Traditional destructive 
readout of the clock state populations in 11) and 12). (b) QND 
probe of the clock transition |1) — |2) using a far detuned 
probe on 11) — |3). The wavelengths are given for Sr. 


achieved for 7 = 1/2 and the Heisenberg scaling of the 
clock precision ^/^/N with the atom number. Such 
scaling for a microwave clock probed on two cyclic 
transitions has been demonstrated in ifTSlISOl . 

The unique energy level structure of alkaline earth 
atoms provides an ideal configuration to implement a 
QND protocol based on a cyclic transition now in an 
optical clock. As a specific example, we consider an 
optical lattice clock operating on the ^Sq (| 1 )) - ^Pq 
(| 2)) transition with N atoms placed in an optical res¬ 
onator. The collective QND readout is performed on 
the cyclic ^Sq (|1)) - ^Pi (|3)) narrow transition (^7.4 
kHz) in Sr using a far detuned probe (Fig. [^b)). The 
clock sequence and the details of the QND measure¬ 
ment are outlined in the Supplementary Material. 

Towards our aim of demonstrating Heisenberg scal¬ 
ing with the number of distant clocks in a clock chain 
we consider first a single lattice clock placed in an 
optical resonator (Fig. [^a)) ISTll . In the following 
we assume that the detuning of the atom and cav¬ 
ity resonances is much greater than the vacuum Rabi 
frequency Q, which in turn is much greater than the 
atomic (7) and the cavity (F) linewidths. 

We consider a standing wave cavity with in¬ 
put/output mirror power transmission coefficients 
Ti, T 2 , single pass intracavity losses jC and the detuned 
probe single pass absorption d^- The cavity power 
transmission coefficient on resonance for small 
is 4Ti T 2 / (Ti + T 2 + 2£ + 2dA) ^ = 4 Ti T 2 / (Ti + T 2 + 
2£)2 [1 - 2FdA/TT], where F = 2tt/{Ti+T2 + 2C) is 
the cavity finesse. Thus, d^, as well as the correspond¬ 
ing phase shift, is enhanced by a factor 2 F/ 7 r. Depend¬ 
ing on details of the experimental realization, the opti¬ 
mal measurement is achieved either in reflection from a 
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single-ended over-coupled cavity with Ti ^ T 2 , 
or with a symmetric cavity in transmission. The atomic 
absorption, d/\ = N/2n at the optimal r] = 1/2, can 
be reduced by using larger A and photon number n. 

Eq .|^is modified in the presence of the cavity. With 
n photons detected during the interaction time the ob¬ 
served probe phase shift consists of two terms: 

(p = + \/ G ! A2F I'kXa- (4) 

The first term is the shot noise of detected light, and 
the second term represents the cavity enhanced phase 
shift. To derive the cavity-based input-output equation 
we multiply both sides of Eq|^with ^/n, 

pout ^ pin p 

= 1 + Vdne~'^^{'y I A)y^alA2Fl7YXA‘ (5) 

Here, Kcav = V dne~'^ / A) a / A2F / n = 
y^dr]ne~^^2F/7T is the cavity enhanced atom-light in¬ 
teraction constant and r]n corresponds to the detected 
photon number n. The relation between the sponta¬ 
neous emission rate in the cavity and in free space 
is Tfcav = VuF/tt for small C. Eq. |^is then modi¬ 
fied with substitutions rj rjcav = 1/2 and Hopt ^ 

K-cav = with 

the optimal value r] = 1/2. Eor the case of large op¬ 
tical depth and/or finesse, AdrinF‘^ j ^ 1, we ar¬ 
rive at the squeezed spin variance irnin = \f^l^cav — 
e7Tl{2dF) oc {FN)~^ achieved for r^cav = 1/2, valid 
for dF ^ 1. The clock precision is then = 

y/27reA/{aF)N~^. Note that our treatment is lim¬ 
ited to aFfA < 27re, otherwise the Holstein-Primakoff 
approximation breaks down when the size of the anti¬ 
squeezed component becomes comparable with N. 

In a realistic design for cavity QED with Sr 
atoms 13^ . we envision F = 10^, length of 5 cm, 
and N = 1000. The atoms in state |1) are collec¬ 
tively coupled to a single mode of this cavity through 
state 13). The bare cavity mode is dressed by the pres¬ 
ence of the 11)-atoms and the resonance spectrum is 
determined by the atom-cavity detuning (A) and the 
collective vacuum Rabi splitting (Q) that depends on 
the number of atoms in |1) (Eig.|^b)). The estimated 
Q = 500 kHz, which is to be compared with T = 29 
kHz and 7 = 7 kHz. The collective cooperativity fac¬ 
tor 0^/(r7) = dF = 1200, leading to an estimated 
20 dB of metrologically useful spin squeezing. Eor Sr 
with a cyclic optical transition, the actual value of A 



FIG. 2. Entangled clock sequence (a) Cavity QED used for 
QND probe of the clock states, (b) The eigenfrequency spec¬ 
trum for 10^ Sr atoms distributed in an optical lattice inside 
a 5 cm cavity with F = 10^. The single-atom effective vac¬ 
uum Rabi frequency is 16 kHz. (c) A cascaded cavity system 
to entangle multiple independent spin-squeezed clocks. 


does not play a fundamental role, but A ^ Q can be 
useful if large values of n are desired. 

It follows from the above discussion that in the 
case of a lossless optical channel connecting M iden¬ 
tical clocks, a collective QND probe of the whole 
network leads to the precision that is a factor of M 
better than each clock, as opposed to uncorrelated 
clocks for which the precision improves by a factor 
of \/M. Eigure |^c) shows an optical probe field 
passing through a chain of successive optical cavi¬ 
ties followed by a single quantum measurement per¬ 
formed at the output. Such interaction generates a 
collective squeezed state of the entire system of M 
clocks. A channel with finite losses can be accounted 
for by the substitution k = ^/4:dFrje~'^/ir Ki = 
with the probe induced decoher¬ 
ence for the i-th clock rn, where describes the op¬ 
tical channel transmission from the i-th clock to the 
detector (subscript cav omitted for brevity). The noise 
of the measurement is still the shot noise of the de¬ 
tected probe whereas the signals due to the spin projec¬ 
tion from all clocks add up, so that S'/A for the chain 

M M 

is ^ ^ MNFir]ie~‘^^^ /tt. 

i=i i=i 

Maximal collective spin squeezing for the chain can 

be found by optimizing this expression, given the 
clock parameters and the channel transmission prop¬ 
erties. Consider, for example, M clocks connected 
with a channel with equal transmission e~'^ between 
each pair of clocks (r^ = {M — i)r and total chan¬ 
nel transmission is / = With the opti¬ 

mal value r]i = ne^'^ jA)‘^a jAFi jn = 1/2, the col- 
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lective S/N becomes MNFm /'^e ^ e = 

i—\ 

_ 1 ) where the fixed 
value of r]i dictates that the cavity finesse is maximal 
for the last clock in the chain, Fi = 

Assuming a sufficiently dense chain of distant 
clocks (r<^l but Mr^l), we reach the precision for 
the chain {S/N)-^ = {MNFm/ 7Te)-^/^M-^\\nt\(x 
(NM) “ ^ I In 1 1 . The expression in parenthesis is limited 
to <^N because the size of the antisqueezed quadrature 
must be Within this limit we obtain Heisenberg 
scaling of the precision of the chain with both N and M 
for any given channel transmission t. For example, four 
clocks probed by QND measurement through a channel 
with t = = 0.5 (3 dB total losses) provide 

precision improvement of 3.1, and eight clocks in the 
same channel give the improvement by 6 . If the ulti¬ 
mate performance of each clock dictates an upper limit 
on N due to, e.g., atomic interactions Ea, a chain of 
entangled clocks may provide an optimal solution. Dis¬ 
tant clocks in a collective SSS may offer an opportunity 
for testing sensitive relativistic effects 1341 . 

EPR entangled clocks. SSS discussed above allows 
for determination of one of the two quantum projec¬ 
tions of the Bloch vector to better than which is 

sufficient for improved clock precision. However, for a 
pair of suitably designed clocks a more intriguing state 
is possible where both projections are defined better 
than this limit with respect to each other. Clock com¬ 
parison can thus run significantly better than the con¬ 
ventional synchronous mode |33l |35l |36l. Such state 
of two Bloch vectors (spins) is a special case of the 
EPR state with the entanglement condition Var{JyiF 
Jy 2 ) + + Jz 2 ) < 2J |29l. It can be re¬ 

alized when the mean spins of the two ensembles are 
oriented in opposite directions, J = Jxi = —Jx 2 , as 
demonstrated for collective magnetic spins E). For 
optical clocks the requirement of oppositely oriented 
mean spins means that the two clocks should be initial¬ 
ized in two opposite clock states (Fig. |^. 

When the Bloch vector describes a pair of states sep¬ 
arated by an optical transition, the conventional QND 
method of generating an EPR state is not applicable be¬ 
cause it would require a direct measurement of the os¬ 
cillations at an optical frequency. However, as demon¬ 
strated for magnetic spin oscillators 13811 . the EPR state 
can be generated by a common dissipation process pro¬ 
vided by forward scattering of indistinguishable pho¬ 
tons that does not involve any measurement. The inter- 


(5s6s) Clock 1 



(5s6s) Clock 2 



FIG. 3. Transitions driving two clocks into an EPR- 
entangled pair. Clock 1 and clock 2 are driven with four 
phase locked classical fields (solid arrows). Forward scat¬ 
tered quantum fields (dashed arrows) generate entanglement 
corresponding to the two clock Bloch vectors being exactly 
antiparallel despite their individual quantum noise. 


action Hamiltonian between two atomic ensembles and 
light that generates an EPR state of the atomic opera¬ 
tors b\ and 62 is iif oc gia^_^b\ -f z/iaL^i -|- + 

r' 2 Ci\.b 2 Fh.c. The first (last) two terms describe the cre¬ 
ation of photon fields a+,a_ and corresponding cre¬ 
ation/annihilation of the collective atomic excitation bi 
( 62 ). For an optical clock these operators correspond 
to the collective excitation generated in the lower (up¬ 
per) state of clock 1 (2). Entanglement is generated 
under the following conditions 13^ : photons scattered 
from the two clocks into mode a+ are indistinguishable 
(same for mode a_) and /ii = /i 2 , = ^ 2 - 

The challenge of realization of such interaction for 
an optical clock transition (or for any collective exci¬ 
tation scheme realized on an optical transition) is that 
due to the selection rules the above conditions are not 
feasible with a standard Raman transition (four-wave 
mixing). It turns out, however, that these conditions 
can be fulfilled using a six-wave mixing process shown 
for a specific example of Sr optical clock in (Fig. |^. 
The use of two-photon driving fields (blue, red, black 
and green solid arrows) allows to fulfill the condition 
of indistinguishability for photons a+ (a_) emitted by 
the two ensembles by choosing the two-photon detun¬ 
ings Sp,6d tohc the same in both clocks and by phase 
locking of the lasers (solid arrows). The condition 
hi = h 2 ^J^i = ^2 ior scattering amplitudes in the 
Hamiltonian can be met by tuning the one-photon de¬ 
tunings (^ 52 , ^Pi^^P 2 - Similar to SS, the degree of 
entanglement scales with the optical depth thus bene¬ 
fiting from cavity enhancement as well. 

An ideal entangled state of this kind corresponds to 
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the two clock Bloch vectors being exactly antiparallel 
(Fig.|3j. This is to be contrasted with the case of SSS 
where the Bloch vector direction is defined better than 
QPN only in the plane in which the squeezing axis lies. 

The EPR state can be used for secret time sharing 
analogous to the quantum key distribution. The clock 
sequence resembles the standard clock sequence (see 
Fig. 1^ in Supplementary Material) with an important 
inset in step (d). At this step the two clock owners ran¬ 
domly choose either to apply or not to apply the 7r/2 
rotation around x axis. They then publicly exchange 
the choice with respect to the 7r/2 pulses, but not the 
results of the clock interrogation. The procedure is re¬ 
peated several times. In close analogy to the quantum 
key distribution we can use the measurements in which 
we have made the same choice of rotations for the rel¬ 
ative time measurements with high precision. Each of 
the clock owners acting separately will achieve a much 
worse precision, compared to QPN, since one half of 
the EPR state is a noisy thermal state. If an eavesdrop¬ 
ping attempt is made, the combined two-clock preci¬ 
sion will be compromised as well. 

Another attractive feature of the EPR entangled 
clocks is the improved capability to check any clock 
disagreement quickly, enabling an efficient approach 
for characterization of systematic effects of an un¬ 
known clock (2) using a well-calibrated clock (1). 

Perhaps the most important application of an EPR 
pair of clocks is for remote sensing. The EPR correla¬ 
tion can be used to map out electromagnetic field from 
site 1 to site 2 remotely. Namely, one can slave clock 
2 to clock 1 by matching the conditions of clock 1 to 
that of 2. In fact, this might be the best tool to explore 
gravitation potential-induced decoherence such as de¬ 
scribed in Ref. 13^ , and it can also serve a potentially 
important role for a future long-baseline atom intefer- 
ometer for gravitational wave detection 1401 . 
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SUPPLEMENTARY MATERIAL 

Quantum nondemolition measurement in an optical 
clock. We consider an optical lattice clock operating 
on the ^So (|1)) - ^Po (|2)) transition with N atoms 
placed in an optical resonator. In a normal clock opera¬ 
tion, the population of the two clock states is measured 
destructively by scattering photons with the strong ^Sq 
- ^Pi transition (Fig. [^a)), permitting a QPN-limited 
probe of the clock transition. A phase sensitive probe 
based on this strong transition was implemented to en¬ 
able a less destructive measurement of the state popu¬ 
lation EH. Here we consider a collective readout of 
the cyclic ^Sq (|1)) - ^Pi (|3)) narrow transition (^7.4 
kHz) in Sr using a far detuned probe (Fig.l^b)). 

An ensemble of N clock atoms can be described by 
the collective pseudo-spin vector J of spin- 1/2 parti¬ 
cles. Jz is defined by the population difference AA^, 
such that: Jz = — N 2 ) = AN/2. Atoms are 

prepared in a superposition of the two clock states by a 
7r/2 rotation of J around the ^-axis of the Bloch sphere 
(Fig. I^b)). The population of |1) is measured by a 
probe detuned by A from the cyclic 11) — |3) transition. 
Note that this probe does not cause redistribution of the 
populations between the clock states, hence Jz is con¬ 
served and is a true QND variable |[T8|[29l|3Tl. After a 
TT-pulse is applied to swap the clock states, the popula¬ 
tion of |1) is measured again. Under this operation the 
effect of the imprecision of the 7r/2 pulse and fiuctua- 
tions of N are suppressed. This QND probe will intro¬ 
duce a Stark shift of the clock transition consisting of a 
mean value and a random shift due to the shot noise of 
the probe (quantum back action of the measurement). 
The former can be canceled by choosing the detuning 
—A for the second measurement of the population. We 
note that the precision of the tt rotation should be bet¬ 
ter than This sequence results in creation of an 

SSS shown as an ellipse in Fig.|^c). Squeezing of Jz 
is then converted into squeezing of the coherence be¬ 
tween the clock states through a 7r/2 rotation around 
the x-axis (Fig. Qd)). The atomic spin is then let to 
process, as in a standard Ramsey sequence (Fig.^e)). 
After a certain precession time, a 7r/2 rotation around 
y is applied (Fig. Qf)), where the population measure¬ 
ment noise is reduced by squeezing. 

Heisenberg scaling for optical clock precision with 
the number of atoms. A coherent spin state (CSS) 
is a product state of individual uncorrelated atoms ori¬ 
ented in the same direction, for example, Jx = N/2, 
with |T^) = + |2)i). The other two pro- 



FIG. 4. Entangled clock sequence. The sequence of oper¬ 
ations of the clock including generation of a spin squeezed 
entangled state. Details in the text. 


jections of J have minimal equal variances allowed 
by the Heisenberg uncertainty relation: Var{Jz) = 
Var{Jy) = Jxl‘^ = A^/4. These fluctuations, re¬ 
ferred to as QPN, and shown as a circle in Fig. [^b), 
pose a fundamental limit to the precision of the clock 
operating on N uncorrelated atoms (271. Introducing 
quantum correlations between atoms allows to reduce 
Var{Jz) to below the QPN limit. For this SSS, under 
the condition that 

Var{J,) < 

^ = VariXA)j < 1 , ( 6 ) 

the atoms become entangled (28l, and the correspond¬ 
ing signal-to-noise ratio for spectroscopy is improved 
by the inverse of ^ 03, which is the spin squeezing 
parameter for metrology. Here we introduce the canon¬ 
ical operators, Xa = Jz/Pa = Jy/ obeying 
the commutation relation [Xa^Pa] = 'i where J is the 
length of the mean pseudo-spin vector, which is the 
atomic coherence. The quantum noise limited clock 
precision defined as the minimal detectable angle of 
spin rotation for the clock sequence is y/Va/r{J^ IJ — 

VW- 

Eq. shows that ^ is determined by the variance 
of the squeezed component of the spin and by the 
mean spin J. The input-output relation for light, with 
similarly defined operators of and P^, assuming 
the QND interaction H oc XaXl (29l is = 

+ tiXA^ To within a factor of unity the interaction 
constant n = ^/dr]e~^, and d is the resonant single 
pass optical depth. Furthermore, y = ndrillA’^ajA 
is a parameter describing spontaneous emission caused 
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by the probe, which leads to the reduction of coherence 
as J = e~'^Nl2. cr is the resonant dipole cross section, 
A is the beam cross section, 7, A the natural linewidth 
and detuning of the optical transition, Udr is the photon 
number in the QND probe. Taking the reduction of J 
into account, we get the SSS with ^ The 


minimal value irnin = — ^e/d oc N~^ is 

achieved for 7 = 1 / 2 , and the optimal tiopt = \/dl2 
valid for d >> 1. The precision of the clock ^/^/N 
then follows the Heisenberg scaling 1/N. 



